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Abstract 

We test a new approach to bosonization in relativistic field theories and 
many-body systems, based on the use of fermionic composites as integration 
variables in the Berezin integral defining the partition function of the system. 
The method appears promising since at zeroth order it correctly describes 
the propagators of the composites, which can be evaluated in a number of 
significant cases. Still to be established is a general procedure for deriving the 
free action of the composites starting from the one of the constituents. To shed 
light on this problem and to explore further features of the method we study 
a simplified version of the BCS model. In this simple case the action of the 
composites can indeed be obtained: whether this result can be generalized it 
remains however to be seen. Yet an interesting property of the wave operators 
appearing in the free actions of bilinear composites already emerges from the 
simple problem we have treated: here the wave operators do not describe 
any time evolution, even though they generate the right propagators. This 
outcome relates to the basic properties of the integrals over the even elements 
of a Grassmann algebra where the composites live, which entails that the 
propagators are no longer the inverse of the wave operators. 



*This work is carried out in the framework of the European Community Research Program 
"Gauge theories, applied supersymmetry and quantum gravity" with a financial contribution under 
contract SC1-CT92-0789 . 
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1 Introduction 



In the present paper we discuss an approach to deal with fermionic composites in 
the framework of the field theory of elementary particles and of many-body systems. 
This approach is based on a change of variables in the Berezin integrals defining the 
partition function of the fermion system, whereby one introduces as integration 
variables the composites one is interested in [0. Although the method is fairly 
general, being suitable for describing trilinear as well as bilinear composites, in this 
paper we confine ourselves to the latter case. 

In the physics of elementary particles fermionic composites appear in general as 
bound states, except when the vacuum is nontrivial due to fermion condensation. 

In the domain of many-body systems bilinear composites are instead related 
to collective excitations which behave like bosons. A paradigmatic example is of 
course provided by the Cooper pairs in metals |2| and in atomic nuclei 0, but other 
well known phenomena of bosonization are met, e.g., in the density fluctuations of 
the electron gas f4[ (described by the Tomonaga model) and in the spin waves in 
ferro-antiferromagnetic metals 0. 

As it is well known, the standard bosonization approach in the theory of many- 
fermion systems is based on the so called Stratonovitch-Hubbard [Q representation 
of the partition function, where the fermion action, quartic in the fermion variables, 
is transformed, via the introduction of an auxiliary bosonic field, into an action 
quadratic in the fermion variables, which allows an exact integration over the latter 
to be performed. If the resulting fermionic determinant can be expanded in a power 
series, the kinetic terms needed to transform the auxiliary field into a dynamical one 
can be generated. The latter field describes then the collective excitations associated 
to the composite. A condition necessary for expanding the fermionic determinant 
is that these excitations should be highly collective; however in general this is not 
sufficient, as we shall see in a specific example. Moreover the method does not 
account for the Pauli principle, which can be a serious shortcoming in atomic nuclei. 
Apart from these limitations, the Stratonovitch-Hubbard transformation is rather 
successful in the theory of many-body systems. 

In the field theory of elementary particles, where the fermionic terms are usually 
quadratic (quartic fermionic interactions being not renormalizable in four dimen- 
sions), the Stratonovitch-Hubbard transformation is of no use. In this context other 
methods have been proposed [0], but they appear to fail in providing a general 
solution to the problem ||: in particular none of them satisfactorily deals with 
composites of low collectivity. 

The approach discussed in the present paper proceeds through the following 
steps: 

i) the Grassmann variables associated to the composites are introduced as new 
integration variables, 

ii) the free action for the composites is next determined, 
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iii) the composite propagators are then evaluated 
and, finally, 

iv) the correlation functions in presence of an interaction are obtained. 

This approach rearranges the usual perturbative expansion in a hopefully ad- 
vantageous way, because it describes the propagation of the composites already at 
the zeroth order. For illustrative purposes we will apply it to a simplified version 
of the BCS theory, the so-called pairing model. The simplest version of the latter 
considers a system of fermions living in an energy level of given angular momentum 
j and interacting pairwise only when the partners are coupled to zero total angular 
momentum. Such a model is a schematic representation of the physics of atomic 
nuclei far from closed shells, where collective excitations associated with a spin zero 
composite occur. Its spectrum is well known from the hamiltonian formalism ||: 
we derive it as well in the path integral formalism, which, as far as we know, is a 
new result. 

We now assess the status of the above program, outlining at the same time the 
content of our paper. 

The first step has been accomplished in refs.[[y], whose results are summarized 
in Sec.0 for the ease of the reader. 

The second step has not yet been fully completed. Indeed first one has to provide 
an unambiguous definition of the free action for the composites. Although the 
natural definition, when it exists, would appear to be the one yielding the free 
propagators appropriate to the spin of the composite, the calculation of the latters is 
far from trivial. For example for a scalar or spinor composite, in order to achieve this 
goal one has to evaluate the correlation function of bilinear or trilinear composites 
with an action which involves the forth or sixth power of the fermionic variables. It 
is a notable feature of the present approach that it allows the evaluation of the free 
propagators in a number of relevant cases, which will be later described. Secondly 
a method has to be devised in order to build the action of the composites out of the 
one of the constituents. 

Clearly in relativistic field theories an action with the above outlined features 
can be added by hands to the fermion one. Indeed in space-time dimensions greater 
than one it is a non-renormalizable, irrelevant operator, which modifies the renor- 
malization group parameters by a finite amount. But one should notice that, first, 
this possibility does not occur in the framework of many-body systems, and, sec- 
ondly, it would be clearly preferable to derive the action of the composites from the 
original fermion one. In Secs.^ and ^ we shall present indeed one way to achieve 
this scope for the simplified BCS model. 

Let us now discuss the correlation functions. We already know [1] that for 
trilinear, spin | composites, the free action is, under certain conditions, the Dirac 
one, but the free action for scalars is still unknown. A candidate, defined on a lattice, 
has been considered in ref. [p]]: its propagator is related, in the continuum limit, to 
the self-avoiding random walk and lends itself, in a certain region of parameters, to 
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analytic approximations || . Rather surprisingly the wave operator appearing in this 
action does not correspond to a true time evolution. Yet a meaningful propagator 
can still be obtained because, unlike in the case of fermions and true bosons, for 
even composites the propagator is not the inverse of the wave operator, as we shall 
show in Sec.||. 

In Secs.||] and |] of the present paper we derive a free action for composites, whose 
wave operator again does not describe any time evolution, although the associated 
propagators are the correct ones. In this connection it is interesting to observe that 
the action of the composites occurring in the pairing model has the same structure 
of the strong coupling mesonic action of QCD @, |l^]. This suggests that the wave 
operator appearing in the latter case fails to describe the time evolution not because 
of the nature of the strong coupling, intrinsically far from the continuum, but rather 
because of a a general property of even composites. 

In Sec.[| we formulate as well the path integral for the simplified BCS model 
in terms of the even variables associated to fermion pairs with zero total angular 
momentum and we deduce the related spectrum. In Sec.| we show how a two-body 
interaction can be handled in a perturbative scheme. For simplicity we confine 
ourselves to the monopole term of the multipole expansion of a general two-body 
force. 

Finally in Sec.^| we present our conclusions emphasizing which parts of the 
present application can be relevant to relativistic field theories. 

2 Change of variables in a Berezin integral 

Let us consider the composites 

£ ^r 2 A n A l2 (i) 

where the A are the generating elements of a Grassmann algebra. The variables 
are characterized by the index of nilpotency n*($/), which is the smallest integer 
such that 

= for n > n*($ 7 ) . (2) 

Our aim is to introduce the ^-variables as integration variables in the Berezin 
integrals defining the composite correlation functions 

< >=Y I l^*d\W h ...$ In e- s ^> x \ (3) 

In other words we must define an integral over the ^-variables in such a way 
that, for an arbitrary function /, 
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Eqs.(|TD obviously cannot be inverted, namely the A cannot be expressed in terms 
of the $. Nevertheless the latter can be considered as new variables in the sense we 
specify below. 

Indeed we first recall that if we have 2Q constituent variables only one func- 
tion (modulo a numerical factor) of the latter exists having a nonvanishing Berezin 
integral. This function is the product of all the A, namely 

A = AiA 2 ...A 2n (5) 

and its Berezin integral is 

J [dX]A = 1 , with [dX] = dX 2 ndX 2 n-i...dX 1 , (6) 

all other integrals identically vanishing. To find out the rules of integration over the 
variables $, we must therefore determine all the functions of the latter which, when 
expressed in terms of the constituent variables, are proportional to A with a nonzero 
coefficient. We call these functions relevant in fact only when relevant functions 
can be built out of the given composites can the $ be introduced as new variables 
of integration. 

Now, since the most general function of nilpotent variables is a polynomial, 
it is sufficient to determine all the relevant monomials, namely the monomials of 
maximum degree. These are referred to as fundamental if they can be expressed 
as products of powers of the $ having the unity as a coefficient and such that 

B n = = wnA, w H ^0. (7) 

The components nj of the vector index n, which characterizes these monomials 
together with the weight Wft, satisfy the constraint 

5>j = fi. (8) 
i 

Notably, any integral over the variables $ can be expressed in terms of the 
weights. In fact a function of the $ can always be expanded as follows 

/(^O = ^ /n®n + irrelevant terms . (9) 

n 

If / is expressed, via the definition of the $, in terms of the A, then its Berezin 
integral is 

J[dX}fmX)]=J2fnWn. (10) 
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This result holds as well when we integrate over the variables providing the 
rule of integration 

J[d$}e H = w H (ii) 

holds for the fundamental monomials, all the other integrals being zero. Note that, 
although in general different expansions 

/(*) = Y + irr - terms = Y fn®n + irr. terms (12) 

n n 

exist, the above equation implies 

J2f™ w * = J2fn w ri, ( 13 ) 

n n 

since both the l.h.s. and the r.h.s. are equal to the coefficient of A in the expansion 
of /(<&) in terms of the generating elements. Hence the value of the integral does 
not depend upon the particular expression for /. 

It should now be clear in which sense relations like (|l|) can be viewed as a change 
of variables, even though they cannot be inverted. In fact in any nonvanishing 
Berezin integral the integrand must be proportional to A, which can always be 
replaced by relevant functions of the composite variables. 

This integration has some distinctive features: 

1) the number of even variables needs not to be equal to half the number of the 
fermionic ones ||. 

2) The integral of a fundamental monomial is unaffected by an arbitrary shift 
$ — >• $ + a. This property, however, does not apply to the integral of an arbitrary 
function: in fact the expansion of the latter in fundamental monomials changes 
as a consequence of the shift. Unlike the Berezin integral, therefore, the integral 
over even elements of a Grassmann algebra is not translationally invariant. Yet an 
equation of motion for the even Grassmann fields can still be derived. 

3) The integral over even Grassmann variables is not invariant under a trans- 
formation from even to even variables such to alter the index of nilpotency. This 
includes Fourier transformations. For this reason the propagator of even Grassmann 
variables cannot be evaluated in the standard way. 

4) The integral of the exponential of a quadratic form of variables of index 1 is 
straightforward, namely 

J [d$*d<l>]e£^** lMl ' J * J = perM , (14) 

where perM is the permanent of the matrix M. This is the reason why the prop- 
agators for even Grassmann variables do not coincide with the inverse of the wave 
operators Mj t j and are difficult to calculate. In fact the permanent of a matrix 
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cannot be evaluated by diagonalization, in turn reflecting the non-invariance of the 
integral under unitary transformations. 

In the present paper the fermionic variables X m are associated to nucleons of given 
angular momentum j and ^-component m. We first introduce then the even variables 
m associated to nucleon pairs with zero component of the angular momentum on 
the spin quantization axis 

<p m = A_ m A m , m > . (15) 

Obviously they have index of nilpotency 1 and only one fundamental monomial 
of weight 1, viz. 

e = <hi2..4i , (is) 

can be built out of them. 

Next we consider the variable $ expressed as a superposition of nucleon pairs of 
zero total angular momentum according to [] 

*= E(-!r m 0-- (it) 

m>0 

The index of nilpotency of <£> is Q and only one fundamental monomial, namely 

= $ n , (18) 

is associated with $: its weight is Q\. 

To perform integrals over the even Grassmann variables when powers of the 
variables $ are present the following identities are helpful: 



J[d<p\<$> n = 5(n - tt)tt\ = J d$<f> n 
J [d<p\<p m <$> n = S(n + 1 - - 1)! = i J 
J [#]0 ni ...0„ fc <P = 5(n + k-Q)(Q- k)\ = - {Q ]_ k + 1) J d^ n+k .(19) 

In fact from the above one sees that w.r. to the integration over $, a product of 
behaves as 



$ fe . (20) 



lr The coupling of the two fermions to zero angular momentum would actually entail a factor 
(— l) J_m /\/2_7 + 1. In this case the pairing interaction appearing in Eq.(39) would be multiplied 
by a factor fi. We have assumed the definition (0) in order to simplify the notations. 
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3 The free action and the free propagator of even 
fields of index 1 

In this Section we introduce the free action for the 0-variables and evaluate their 
propagator. The partition function of a free system of fermions reads 



= J[d\*d\]e~ So (21) 



where 

_ T 

[d\*d\] = J] II d\ m (tyd\ m {t) (22) 



m t=l 



is the volume element and 



So = rJ2J2 **m(t)(Vt + e m )X m (t - 1) (23) 

t=l m 

is the euclidean action. In the above r is the time spacing, Vj is the discrete time 
derivative 

V,/(f) = ^[/(* +!)"/(*)] (24) 

and e m is the fermion energy. Note that the variables at the time t — are related 
to the ones at t = T by antiperiodic boundary conditions 

A m (0) = -A m (T) . (25) 

Let us now introduce the even elements of the Grassmann algebra, which describe 
pairs of fermions with J z = 

* m (t) = \- m (t)\m(t), = ^m-Jt). (26) 

We want to write the partition function as an integral over these variables 



= J [d(f)*d(f)}e- So . (27) 



As a first step to express in terms of the 0's let us write the free action as 
follows 



So = e E {^(*)A m (o + y m m- m (t) 

t=l m>0 

- x m [\* m (t)\ rn (t - 1) + A*_ m (0A_ m (* - 1)] } (28) 



8 



where 



■Em 1 T~e m • (29) 

Next we express the exponential of (minus) the free action as the product of 
exponentials (associated with different times and different absolute values of m) 
and expand each exponential accounting for the nilpotency of the variables. We get 



T 



So 



n n { I 1 + cw<u*) - (A^(t)A m (t) + \*_ m (t)\„ m (t)) 

t=l m>0 



(30) 



[l + x 2 m cj> m {t)cj> m {t - 1) + x m (X m (t)X m (t - 1) + A* m (t)A_ m (£ - 1))] } 
or, by explicitly performing the products, 



-So 



= n n {i + + xi^t^ut - 1) 

t=l m>0 

+ x m (j)* m (t) [A m (t)A_ m (t - 1) + X m (t - l)A_ m (t)] 



(31) 



A^(t)A m (t) + A* m (t)A_ m (t)] + a; m [A^(t)A m (t - 1) + A* m (t)A_ m (t - 1)] } 



Now one easily realizes that the terms in the last row of the r.h.s. of (plD contribute 
neither to Zq nor to the expectation values of operators involving only the <p m and not 
the A m variables. Hence they can be neglected in the evaluation of the propagator 
of the 0. 

In order to get the partition function, one has to identify in ( j31~l) the terms which 
generate monomials involving all the <fi* and 0's. One obtains 



{T T 
1] 4> m {t)<l>m{t) + I] X 2 m <P* m {t)^mit - 1) 
t=l t=l 

+ II XmFJt) [A m (t)A_ m (t - 1) + X m (t - l)A_ m (*)]| , (32) 
t=l ) 

where the integrals are straightforward and yield, in the limit T — > oo, the well 
known result 



(33) 



m>0 m>0 

Let us now evaluate the correlation function between two composite variables 



< 0„(t 2 )0;(ti) >= ir I [d<P*d4>]<p n (t2)<P* n (ti)e~ s ". 



(34) 



Proceeding as before we find 



1 p -(/3 2 -/3i)2e„ 

z o m ^ n [1 + e-^j 



for ti < i 2 and 



1 p -(/3+/3 2 -/3i)2e n 



for ti > t 2 - Note that Eqs.(^) and ( pq ) are just the squares of the corresponding 
correlation functions for the A variables. However, here the existence of a back- 
ward propagation is an artifact due to antiperiodic boundary conditions. We must 
therefore take (3 3> |/3 2 — A I m order to get rid of this artifact. 

Note also that now the term linear in x m in Eq . fl3T|) does not contribute to the 
numerator. Accordingly it will not affect the energy of the system and we can 
assume as free action 

# =-££ [4>*m(t)<f>m(t) + X 2 m4>*m(t)4> m (t - 1)] . (37) 
t=l m >0 

This action is quite different from the true boson and fermion ones, because in 
the continuum limit it does not contain any time derivative. As explained in the 
Introduction we conjecture that this is a general feature of the free actions of even 
Grassmann fields. 

Note that if the partition function is calculated with fl37|), one gets a wrong result 



na+^ T ) — n [i+e- 2/3 -] . (38) 



m>0 m>0 

We can therefore use to evaluate the spectrum and the ratios of correlation 
functions, but not their absolute normalizations. 

Finally we should emphasize that, although the variables <ft are introduced to 
bosonize the system, they are completely equivalent to the A and can be used there- 
fore to evaluate also the propagator of a single nucleon, as shown in the Appendix. 



4 The pairing interaction 



In this Section we consider a system of fermions interacting via a pairing potential, 
whose action reads 



S 



T ( 

t=l lm>0 



(39) 
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where 



$(*) = ]T(-ir m A_ m (t)A m (t) (40) 

m>0 

describes a pair of fermions with zero total angular momentum. We have included 
only the term sjf* from the free fermionic action, because the other terms cannot 
contribute to correlation functions of the $. Moreover we deal for simplicity only 
with the degenerate case, i.e. e m = e, x m = x. 

The spectrum of the system described by the action (|3~9"D , restricted to states of 
n spin zero pairs of nucleons (the so-called seniority zero states), is |9j 

E n = n[2e-(n-n + l)g p ] (41) 

and it is characterized by the parameter 

e-ffi+iJat. (42, 

We will confine ourselves to the case £ < 1, since then the ground state is triv- 
ial. Worth mentioning is that, even in presence of high collectivity, insured by 
a large value of Q, the fermionic determinant arising after the application of the 
Stratonovitch-Hubbard transformation can be expanded in a power series [§] only if 
£ <C 1. Our method is not affected by such a restriction. 

The correlation functions of n spin zero pairs of nucleons are 

< $ n (t 2 ) [$*(t0] n >=\\ Wd4>W{t 2 ) [$*(t!)] n e- 5 . (43) 
Their evaluation is simplified by introducing the variables 

B(t) = j2Mty[ ( p m (t) + x 2 (i )m (t-i 

m>0 

= £G»M, (44) 

m>0 

in terms of which 

S = J2[-B(t)-g p r^m(t-l)) 

t=i 

T 



E S(t) (45) 



t=i 



and 



< $«(t 2 ) [$*(t0] n >= ^ / [dW]* B (t 2 ) [**(«i)] B n E (46) 



f=l r =0 



/[#*#]^ 2 )[$*(ti)] n nE3E flV r - fc (*)(^) fe [^(0^-i)] fe . 

J t=l r=0 r ' fc=0 V / 
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The powers of B appearing in (|46f) are given by 



B r {t) 



E C m(t) 

m>0 



E 



C^(t)...C k n -(t) 



fci+...fcn=r 

r! E ^ (*)-*(*). *i = 0,l 

fei+...fen=r 



(47) 



where the restriction on the indices stems from the nilpotency of the variables C. 
Now each of the ( n j terms (the number of combinations of class r of objects) 



in the r.h.s. of Eq.(47) provides the same contribution. Accordingly in the integral 
B r behaves as 



~r!( n W)...#(t)£( r V^^ 



s=0 



, s , 



due to the symmetry w.r. to permutations of the 0's. Because of this symmetry, as 
far as the integration over the $ is concerned, the rule (|20|) holds, so that 



B r (t) ~ Jj^ [**(*)]' E Q (« - *)! [« - (r - s))\x^ntW- s (t - 1) . (49) 



Now, in relation to the integration over $*, the correlation function is non- vanishing 
Ld only if a power [<3>*(t)] n 
will not vanish only when 



if and only if a power [$*(t)] n appears in the integrand for each time t. Therefore 



r = Q for t ^ t\ 
r = tt-n for t = t\. (50) 

On the other hand, in the integration over the variable $ at the time t = t 2 , in order 
to recover the power $ n (t 2 ), one should single out from B Q ~ k (t 2 ) a factor $ n_n (i 2 ). 
It is then an easy matter to obtain 

J d$*(t 2 )d^t 2 )$ n (t 2 ) [ ~ S ^ ] = a$ n (t 2 - 1) , (51) 

where a is a number which will be evaluated below. This shows that the integral 
over & n (t 2 ) is proportional to $ n (t 2 — 1), so that, by performing the integrals in the 
range t\ < t < t 2 we get 
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* 2 ^S(t)f 

In the limit of vanishing r one has 



J[d^(t)dm] tl<t < t2 ^ n (t 2 ) n [ -^r- = op-^fa) . (52) 



lim a t2 ~ h = lim oft*-*^ = lim e &-*i)/ T * a = e ~E{h-h) 
where the physical time t — rt has been introduced and 

£ = lim f — In a) (54) 



has the significance of the energy. The coefficient a needs to be calculated to first 
order in r only. In turn this entails 

[S(t)] n 



m = ^Eg)B°-*(*)^)*[*-(t)*(t-i)]* 

~o ^[^W + n^B - 1 ^*)^-!)]- (55) 



By inserting fl4"9] ) into (|55p one then gets 
[-S(t)] n [$*(t)] Q 



2 



{$ sz (t) (56) 



+ [x 2(n ~ s) + (s + - s)&,t] $ s (t)$^ s (t - 1)} 

and since 

x 2(n " s) ~ o 1 - 2(fi - s)er (57) 
one can finally cast (|56|) in the form 



+ J2 [1 - (2(0 - s)e - (s + l)(fl - % p ) r] *'(t)* n -(t - 1)} . 

s=0 



Therefore for the relevant integral appearing in the correlation function one obtains 



J d^(t 2 )d^(t 2 )^ n (t 2 ) [ ~ o $"(t 2 - 1) {1 - n [2e - (fi - n + 1)^] r} , 

(59) 

from where it follows 
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a = l-[2e-(Q-n + l)g p ]rn (60) 

and 

< $ n (t 2 ) [$*(tl)] n > ~ e -(ft-A)[2e-(n-n+l) ffp ]n (g^ 

for the correlation function. Hence for the energy spectrum of the system one 
recovers the well known formula (41). 

From the above calculation the free action for the field $ is easily deduced. In 
fact, setting g p =0 in Eq.(^), we see that this action can be written as 



sf> = -s> £>*wr£ r ( n -') | ^( f -')i | jM ff « r . (( - 1 ) 

t=l lr=0 r - s=0 W l iZ V J 

(62) 

Clearly the free correlation functions, which turn out to be independent of the 
total angular momentum of the pair, are then 

< $ n (t 2 ) [$*(tl)] n > r ~ e -(&-ft)2en _ ( g 3) 

We have evaluated the logarithm for a few values of Q, always obtaining an 
action of the form 

5! $) = -EE« s rWr[$ s (t)+^ s (i-l)] , (64) 

t=l s=l 

a s being a numerical coefficient. It is indeed remarkable that Sq^ displays the 
same structure of the mesonic action appearing in the strong coupling limit of 
QCD [|7|, llPf . Such a structure is characterized by the absence of terms of the type 
[$*(t)] n [$(t)]"~ s [$(t — l)] s (s ^ 0), due to the occurrence of nontrivial cancellations. 
Obviously the total action can be recast in the form 

SW = sP-g p Tj2$*(t)m-l) . (65) 
t=i 



5 Perturbation theory 

In this Section we show how interactions can be handled in a perturbative scheme 
considering, as an example, the monopole term of the multipole expansion of a 
quartic interaction. Here the quartic character of the interaction is irrelevant, a 
similar calculation could be performed for a quadratic one, and higher multipolarities 
can be dealt with in a similar fashion. The associated action can be written as follows 
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tR t 

t=l m^m' 

where R is the radial integral of the monopole term of the interaction, and can be 
trivially split into two pieces: one expressed only in terms of the <fi variables and the 
other in terms of the odd Grassmann variables A according to 



Sj = St + S 2 , (67) 
S 1= rfl^CWa-l) (68) 

t=l m>0 

tR t 

5 2 = -yEE E \Ut)\* m ,mm>(t-l)\ m (t-l) . (69) 

t=l rn m'^rn 

For the partition function of the fermion system we shall accordingly have 

Z= f[d\*dX]e- s (70) 



where 



and 



with 

S = So + Sl (71) 

Now in general the terms in the third row of Eq . fl3lD , which expresses So, can 
no longer be neglected because the (^-variables could be reconstructed utilizing terms 
arising from the expansion wr to Sj. 

Thus, in order to be able to exponentiate the whole expression appearing in 
Eq. (|31~D , we introduce the even variable e m (t) of index 1 and rewrite the exponential 
of (minus) the action of the free fermions in the form 



where 



A So = E e m(t) {\* m (t) [Am(t) - x\ m (t ~ 1)] + X*_ m (i) [X- m (t) - xX- m (t - 1)]} 

t,m>0 

-x<f>* m {t) [A m (i)A_ m (f - 1) + X m {t - l)A_ ro (t)] , (73) 

setting e = 1 at the end of the calculations. ASo should now be treated as a 
perturbation together with Si. In each term of the perturbative expansion the A's 
must be combined into the 0's, the integration over the must be performed and, 
finally, the e must be set equal to 1. It remains to be seen whether the expansion 
wr to ASo is permissible in the absence of a small parameter. 
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In the example we are considering we have, in first order of perturbation theory, 



e- s ~e- s o* ) e- AS °(l-S 1 -S 2 ) 



(74) 



The related partition function is thus 



where 



,(e) 



[l + X 



r(e) 



2T\fl 



Zi 



1 



-2/3e 



(75) 



(76) 



and the correlation function of the <p m , which turns out to be independent of m, 
reads 



< ^(tO0 m (t 2 ) >= 



(77) 



Go (ti,tz)/Z being the free pair propagator given by Eqs. (|35|) and (|36|). 

Now we should first notice that odd powers of ASo never come into play: indeed 
the related Berezin integral vanishes owing to the impossibility of reconstructing a 
complete product of the <ft variables. It can moreover be checked that, for the specific 
interaction we are considering, ( f7S|) and ( [77]) are not contributed to by powers of ASo 
larger than two owing to the nilpotency of the variables and e. Thus, in addition 
to Si, only the term quadratic in ASo has to be considered, and this justifies the 
expansion of the exponential wr to ASo- 

Concerning Si one gets 



Z 1 = -tTRQ(1 +x 2T ) n ~ 1 x 2 ^ 
for the partition function and 

G?{t u t 2 ) = -rRx 2 ^~^x~ 2 {l + x 2T ) n - 1 



(t 2 -tO+T(j2-i; 



x 



2T 



l + X 



2T 



(79) 



for the correlation function when t\ < t 2 - 

More delicate is the analysis of the last terms in (|75| ) and ( [77] ) , which brings into 
play ASo via the coupling to S 2 - For these it is indeed possible to reconstruct the 
complete chain of 0's; one thus obtains 



Z 2 = -2rTRQ(Q - 1)(1 + x 2T ) n - 2 x 2 ( 2T -V 



(80) 



and 



gWfah) = -4rR(Q-l)x 2 ^-^x 2(T ~ 1 \l+, 



„2T 



«-<0 4< !i - 2 'TT^ 
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However, in the above expressions the terms proportional to x 2T ~ e~ 2l3e provide 
a negligible contribution, according to the discussion following Eq. (p6|) . As a conse- 
quence one simply gets for the pair correlation function in presence of the interaction 
(|66|) and for t\ < t 2 the expression 

< C(*i)<U* 2 ) I + e -2, e > (82) 

which corresponds to a constant shift of the pair energy 2e — > 2e + R. 

If we now consider the correlation function of two pairs, then the term S 2 does 
contribute through its coupling to AS , yielding: 

< <PUtlWn{tl)<Prn{t 2 )<P n (t 2 ) >= 1 [g^ (t h t 2 ) + g[ 2) (t la t 2 ) + (h , t 2 )] . (83) 

In the above, if t\ < t 2 and again neglecting the contributions proportional to x 2T , 

gW(t 1 ,t 2 ) = {l + a* r ) a - 2 x 4 to-V (84) 
g[ 2) (h,t 2 ) = -2Rr(t 2 - + x 2T f- 2 x^-^x- 2 (85) 

and 

g?\t 1 ,t 2 ) = 2g[ 2 \t 1 ,t 2 ) . (86) 

By taking the limit r — > one finally obtains 

e -(/3 2 -/3i)2(2e+3 J R) 

< 4>* m (tl)K(tl)4>m(t2)cf>n(t2) >^ ^ + e _ 2/3e)2 , (87) 

which amounts to the shift e — > e + |i? for the single particle energy with respect 
to the non-interacting case. 

From the above results the perturbative corrections to the spectrum of a system 
of one and two spin-zero pairs follow. Indeed we have for the correlation function 
of one spin- zero pair 

< ^(hWh) >= £ (-l) 2 ^ m ™< 0^(^)0^^) >, (88) 

mi,m,2>0 

which, the correlator of the <p m being diagonal in m and m-independent, entails 

<$*(t 1 )$(t 2 )>=fi<0^(t 1 )0 m (t 2 )> . (89) 
By a simple extension the correlator 

< [$*(ti)] 2 [$(t 2 )f >= 2Q(Q - 1) < C(ti)0;(ti)0 m (t 2 )0„(t 2 ) > (90) 
is also obtained. 

Clearly this scheme can be easily generalized to obtain the correlation functions 
when both the monopole and the pairing interaction are acting. 
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6 Conclusions 



In this paper we have tested a new approach to bosonization to be used in the 
framework of the theory of relativistic fields and of many-body systems. A new ap- 
proach seems to be required especially in the context of quantum field theory, where 
the Stratonovitch-Hubbard method does not help. The present scheme, hopefully, 
might indeed offer some advantages because first it describes correctly, already at 
the zeroth order, the propagation of composites and furthermore because the intro- 
duction of the composites as integration variables is expected, on general ground, 
to be helpful in fastening the convergence of the perturbative expansion through a 
rearrangement of the latter. 

The formalism utilized in the present work has been developed in previous papers 
where the free propagators have also been evaluated in a number of other relevant 
cases. The problem still waiting for a general solution is how to derive the free action 
of the composites starting from the one of the constituents. In order to explore this 
point and to study additional relevant features of the method, we have studied a 
simplified version of the BCS model. In this simple context it has been possible to 
fully work out the new bosonization scheme, obtaining the spectrum, already known 
in the hamiltonian framework, in the path integral formalism. 

In particular we have been able to derive the exact free action for the composite 
starting from the one of the constituents. This is the first example, to our knowledge, 
where such a derivation has been performed without recurring to approximations. 
It is presently unclear whether our procedure can be generalized to relativistic field 
theories or to many-body systems of greater complexity than the degenerate pairing 
model. We however believe that the structure of our free action is not restricted 
to the spin-zero composite of the pairing model, but is likely to be shared by the 
actions of all even composites, in particular as far as the wave operator is concerned. 
A support to this conjecture is provided by the form of the action proposed in ref. 
! J J and also by the mesonic action in the strong coupling limit of QCD. 

Finally a comment concerning the role played by the pieces of the free action still 
expressed through the original odd Grassmann variables is in order. Terms of this 
kind are indeed expected to be present in any realistic problem and, as a rule, they 
occur without a small parameter in front. In this connection it is of significance our 
finding in the context of the pairing model: here these terms can only be effective 
when an interaction is acting among the constituents and, accordingly, they can be 
treated as a perturbation together with the latter. Also in this case we are inclined 
to believe this to be a general feature of all the actions expressed through the even 
elements of a Grassmann algebra. 

In conclusion the present study shows that the approach to deal with composite 
fields, originally developed in ref.0, can be successfully applied to nonrelativistic 
many-body problems whenever the introduction of complex building blocks (the 
composites) is hinted by the nature of the system. In particular, and importantly, 
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the method allows the derivation of the free action for the composites and appears 
worth to be explored in the context of more realistic situations. 



Appendix 

In this Appendix we want to show that, since the introduction of the is a change 
of variables, we can use them also to evaluate the propagator of a single nucleon 

< A n (t 2 )A;(t 1 ) >= i- / [rf0*rf0]A n (t 2 )A;(t 1 ) e - 5 ° . (91) 

To this purpose let us start by considering the case t 1 < t 2 and n > 0. We shall 
take into account only the term \m\ = |n| in ( |31| ) since the other terms are canceled 
out by corresponding terms of Zq in the denominator. In order to construct 0*(ti), 
we must combine A*(ii) with the terms containing \*L n (t) in the t = t\ term of the 
third row of (|3TD, which yields 



hA-n(ti) + x n X^,{t x - 1)] . (92) 
The two terms of the above expression are the only insertions arising from the third 



row of (0) which, when combined with the terms in the first row at t ^ t\, give rise 
to two "chains" containing all the <fi and <fi* and thus yielding a non-zero Berezin 
integral. The first chain is (a minus sign comes from the interchange of A* (ti) and 
A„(* 2 )) 



- [0;(i)0„(i)] ... [&(ti - i)Mh - m-PMX-nih)] 

X [X n </>* n (ti + l)A„(ii)A_ n (*i + 1)] -. [X n <j)* n (t 2 )\ n (t 2 - l)A_ n (t 2 )]A n (t 2 ) 

x [<(t 2 + l)0 n (t 2 + l)] ... [<(T)0 n (T)] 

= ILKWMt) , ( 93 ) 

t=l 

while the second one reads 



- {Xn^^X^iX-nfi) ... [x n (j)* n {t\ ~ l)X n {tl - l)A_„(ti - 2)][x n (j)* n {tx)X- n {t 1 - 1)] 
X [Xl4>l{t x + l)Mh)} - (^0;(t 2 )0 n (t 2 - l)]A n (t 2 ) 

x [x n (f>* n (t 2 + l)X n (t 2 + l)A_ n (t 2 )] ... [x n( f>* n (T)X n {T)X- n (T - 1)] 

= x^xf^xi-^ n K(t)ut) = xr ti+T n , (94) 

t=l t=l 
where the antiperiodicity condition (|25| ) has been used. Therefore 
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<\ n (h)K(t 1 )>=-x t rr h (l+x T n ) I] [l + ^] 2 rt TT^Fe • ( 95 ) 
Similarly, for ti > t 2 one gets 

< A n (t 2 )A*(ti) >= l^- tl+t2 (l +^) II [1 + ^ ~ ( ~TT P W~ ■ (96) 

|m|^|n| 
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